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SENSITIVITY OF A GAS FILTER CORRELATION 
INSTRUMENT TO VARIATIONS IN OPTICAL BALANCE 

By 

Harry D. Orr, III and Shirley A, Campbell 


ABSTRACT 

The sensitivity of a simple Gas Filter Correlation Radiometer (GFCR) to varia- 
tions in the instrument ^ s optical balance parameter, has been studied 

theoretically. 

A comnuter program was used to simulate the response of the GFCR to changing 
pollutant levels of CO, SO2. CH^, and NH3 in two model atmospheres. 
Positive and negative deviations of of magnitudes 0 . 01 , 0 . 1 , and 1 percent 

were imposed upon the simulation and the resulting deviation in inferred 
pollutant concentration were determined. 

For the CO and CH4 channels, and the higher pressure cell of the NH3 
channel, the deviations are less than ll 2 percent for the deviations in 
of lo.l percent, but increase to significantly higher values for larger devia- 
tions. For the lower pressure cell of the NH^ channel and for the SO2 
channel, the deviations in inferred concentration begin to rise sharply between 
0.01 and 0.1 percent deviation in suggesting that a tighter control on 

may be required for this channel. 

I. INTRODUCTION 

The concept of measuring the total burden or column density of trace gases by 
nondispersive infrared techniques is currently an area of much interest. One 
variety of this technique, the gas filter correlation radiometer (GFCR), is 
being explored for both aircraft and satellite applications to the measurement 
of pollutant gas concentrations {refs. 1 - 2 ). The operation of the instrument 
depends not only upon the various parameters describing the underlying surface 
and atmosphere, but also upon those describing the instrument's physical state 



— for example, the temperature, emissivities and reflectivities of various 
optical components. Moreover, the reliability of the data from such an in- 
strument is strongly affected hy the stability of these parameters. This re- 
port presents the results of a theoretical study of the sensitivity of the 
GFCR to one such parameter, the vacuum - cell aperture transmission, T^, which 
is used to mechanically balance the optical paths of the instrument under 
conditions of no pollutant gas. The role of this aperture is described in the 
following section. A list of symbols is given in Appendix B. 

II. THEORY OF OPERATION OF GFCR 

Detailed descriptions of instruments based on the gas filter correlation 
technique have been presented elsewhere (see, for example, refs. 1«2), In this 
work, only the basic elements of a prototype instrument will be described. 
Figure 1 is a schematic diagram of a simple GFCR. Radiation from the earth s 
surface and intervening atmosphere enter through the front window. One half of 
the energy is blocked by a chopper blade; the other half is passed throiigh a 
lens, then through a cell containing a sample of the gas of interest, and is 
focused onto a detector. One half cycle later, the chopper blocks the energy 
on the right hand side and passes it on the left hand side, which contains an 
empty cell and an adjustable aperture. The difference between the two signals 
can be expressed as a differential radiance, M, proportional to an instrument 

voltage, AV ; 

AV a AN a 

where E{v) is the incident radiation at wavenumber v, C(v), and C'(v) are 
the transfer functions of the optical components of the right and left hand 
paths, respectively, Tq('^) is the wavenumber-dependent transmittance of the 
reference gas, and is the wavenumber- independent transmittance of the 

aperture of the vacuum cell. The filter accepts radiation in the bandpass Av 
The signals also may be combined to give an integral’ radiometL^ measurement: 

V a N a I E(v) [C(v) 'Tq(v) + cMv) L^] dv (2) 

’'Av 


/ 

-'Av 


E(v) [c(v) - C»(v) 




dv 


( 1 ) 


2 



on 


When none of the gas of interest is present outside the instrument, the de- 
pendence of E on V is very weakly correlated with the dependence of 
V. In practice this condition is achieved hy alternately placing a hot and 
then a cold hlackbody at the entrance window. The change in AV is then 
minimized hy mechanically adjusting the aperture transmission, This 

procedure is commonly referred to as "halancing” the instrument. When the gas 
of interest is present outside the instrument, in the atmosphere or in a cali- 
bration cell, the dependence of E on V is altered by absorption and emis- 
sion at the same wavenumbers as are active in the instrument gas cell. This 
correlation produces an imbalance in the signal strengths between the two op- 
tical paths and causes the output signal to deviate from its balance value. 
At the same time, the integral radiance signal, V, is strongly related to the 
brightness temperature of the background surface. 

III. DATA REDUCTION PROCEDURE 


Since it is not currently possible to simulate accurately the infrared activity 
of the earth’s atmosphere in the laboratory, theoretical techniques must be 
used to relate the voltage output of the instrument to the gas burden measured. 
The primary tool used in this effort is a line-by-line radiative transfer 
program "POLAYER" (ref. 3). For the purpose of this program, the background 
surface (either the earth or a laboratory source) is described by a temperature 
and emissivityi the intervening atmosphere is specified in the form of one 
or more layers with a thickness, pressure, temperature, and concentration of 
gases. The instrument is described in terms of the temperature and spectral 
transmitt ances of the various components of each of its optical paths. Then, 
the radiation reaching the detector along each of the optical paths is 
computed by integrating over Av the surface radiation, as modified by any 
absorbing spectral lines in each layer, plus contributions by any emitting 
spectral lines in each layer. This procedure is employed in three stages. 

First, the theoretical value of is computed which will produce, equal 

values of AN for the known laboratory sources with no gases between the 
sources and the instrument. Next, using this value of T^, a description cor- 
responding to the calibration cell used with the instrument is inserted into 
the program. The results of this calculation establish a relation between 
the voltages measured by the instrument and the corresponding values of 
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radiance. The real atmosphere with its inherent non-homogeneity, cannot he 
simulated by such cells; thus in the third stage a description of the real 
atmosphere, to whatever accuracy it is known at the time of measurement, is 
inserted. The response of the instrument to the radiances associated with a 
range of potential concentrations of the gas of interest is computed. Figures 
2-5 show the responses of the instrument to the standard atmosphere used in 
this study. To infer the actual concentration, the voltage calibration from 
stage two is used to convert the instrument output to radiance units (e.g., 
watts/cm sr); then the corresponding gas concentration is found by interpola- 
tion on the response curve generated in stage three. 

IV. METHODOLOGY 

In the present work, the sensitivity of the GFCR was studied in four spectral 
regions, corresponding to wavelengths that might be used in measurement of 
carbon monoxide (CO), methane (CHj^), sulfur dioxide (SO^)! und ammonia (WH^), 
for two pressures of the instrument gas cell (one in the case of SO^) • The 
details of these conditions are given in table 1. 

To use the POLAYER program, two atmospheric models were constructed. In each 
model, the region from the surface to 22 km was divided into five layers. The 
layers were chosen so that each contains an approximately equal column density 
of gas (molecules). The details of the layers are given in table 2. 

Model 1 is based on the U.S. Standard Atmosphere, 1962 (ref, The mean 

temperature of each of the five layers was determined by weighting the U.S. 
Standard temperatures with the corresponding density values. The vertical dis- 
tribution of water vapor was taken as that of the Gutnick Model (ref. 5)» 

Model 2 is based on the U.S. Standard Atmosphere Supplements, 1966, 1+5®N. 
latitude for July (ref. 6). The water vapor distribution was talten from the 
midlatitude summer atmospheric profile of AFCRL Optical Properties of the 
Atmosphere (ref. T). 

For both models, a surface emissivity of .98 was used; a night'ciiue measurement 
was assumed, hence no solar contribution was included. For Model 1, a ground 
temperature of 288K was used. While the concentrations of all other species 
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were held constant, the concentration of the pollutant gas of interest, uniformly 
distributed, was varied in ten steps below and above its nominal or background 
level. For Model 2, a groiind temperatiire of 296K was used. The details of 
the distributions are given in table 2 . 

For each instrument channel and atmospheric model combination, seven response 
function curves were computed, one corresponding to the standard or correct 
value of T. and six other curves corresponding to variations in of 

^. 1 /^, and ^. 01 ^ from its standard value. 

To simulate measurements, four values of AN 'were then selected, and the 
corresponding values of pollutant concentration were determined by interpolation 
on each of the seven curves. For each of the four AN values, the differences 
between the indicated concentration in the standard case and the concentrations 
found in the six varied cases were calculated. 

V. RESULTS 

The results of the procedure described above axe plotted in figures 6-13 and 
tabulated in Appendix A. Each figure, corresponding to a particiilar species, 
cell condition, and atmospheric model, shows the percent deviation of the 
inferred concentration as a function of the percent deviation of from its 

standard value. 

For all cases, the relative behaviors of the deviations are similar. As 
is decreased from its standard values, the percent deviations in inferred con- 
centration become increasingly negative; as x^ is increased, the deviations 
become positive. The rate of change of the deviation also appears to increase 
with deviations in x^. For CO, CHj^, and cell 1 of the NH^ channel magnitudes 
of deviations of inferred concentrations \indergo a sharp rise as the deviations 
in X^ exceed ^0.1^. For cell 2, of NH^ and SO^ the rise commences at even 
lower deviations in x^. 

The deviations for CO, CHj^, and NH^ for a given atmospheric model are lower 
for the cell with the higher concentration. The reason for this behavior can 
be seen in figures 2, U, and 5. The response function for the cell with the 
higher concentration has less of a tendency to saturate in the range of 
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atmospheric concentrations considered. Therefore, the inferred concentration 
at any signal level is less sensitive to the small changes in the response 
function induced by variations in 

For a given species and cell concentration, the deviations in inferred concen- 
tration are of smaller magnitude for the Model 2 atmosphere than for the Model 
1 standard. Although not illustrated in the figures, the reason for this 
situation is similar. The Model 2 atmosphere, with higher surface and atmos- 
pheric temperatures than Model 1, produces generally higher signal levels and 
leads to response functions with less of a tendency to saturate than those 
arising from Model 1. These ’’stiffer" response functions lead to lower devia- 
tions in inferred concentration with changes in 

VI . CONCLUSIONS 

From the preceding discussion, it is clear that the sensitivity of the GFCR 
to variations in the balancing aperture is not a simple relationship. Not only 
do the deviations depend in a very nonlinear way on the magnitude of the devia- 
tion in but also upon the particular spectral region or species con- 

sidered, the concentration in the sample cell of the instrument, and the details 
of the surface and atmosphere over which the instrument is operated. In addi- 
tion, for a given deviation in T^, there does not seem to be a consistent 
relationship between the level of inferred concentration and the deviation to 
be expected. In other words, for a particular T^, it is not possible to 
predict whether the greatest deviations in the inferred concentration will 
occur at high, intermediate, or low pollutant levels. 

The complex dependence of the deviations notwithstanding, some general con- 
clusions concerning optical balance from an operational point of view can be 
reached. 

For the CO, CHi^, and the hi^er pressure cell of the NH^ channel, the 
deviations in inferred concentration will be kept within reasonable limits , 
i.e., to 12 percent, if is maintained to within. lo.i percent cf its 

nominal value. In view of the consistently high deviation observed for the 
lower pressure cell of the NH^ channel, consideration of a higher partial 
pressure in this cell seems warranted. It should be noted that the above 
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conelusions apply to the ranges of pollutant concentrations considered in this 
study and may not be valid for ranges significantly different from them. Also, 
it is evident that atmospheric and surface conditions which tend to reduce the 
signal levels; e.g. , lower surface temperatures, and more nearly isothermal 
temperature profiles, will tend to increase the deviation expected at ajoy given 
deviation in For the lower pressure cell of the NH^ channel and for the 

SO 2 channel, unacceptably large deviations are much more likely to occur with 
io.l percent tolerance; on x^. A rough estimate of the limits on the optical 
balance required to give the same probability of not exceeding about ll2 percent 
deviation in concentration would be ^0, 03 percent. The same qualifications, of 
course, apply as did to the other species; in fact, the sensitivity of the 
required tolerance to surface and atmospheric parameters may be even stronger 
in the case of the lower pressure cell of NH<^ and of SO 2 than for the other 
species considered, 
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Center Wavelength (micrometers) 

Bandwridth (micrometers) 

2 

Cell Pressure (N/m ) 

Cell Length (cm) 

Partial Pressure (N/m^) 

Cell 1 
Cell 2 

Cell Temperature (k) 

Balance Temperature (K) 

Cold 

Hot 

Standard Aperture Transmission (t^) 
Cell 1 


Cell 2 


Table 1 


INSTRUMENT PARAMETERS 


CO 

OJ 

0 

to 

»«3 


4.66 

8.55 

11.2 

8.0 

4.50 - 4.32 

8.11 - 8.99 

10.43 - 11.97 

7.53 - 8.1*7 

1.013 ^ 10^ 

1.013 X 10^ 

1.013 X 10^ 

1.013 X 10^ 

1 

1 

1 

1 

k 

3.546 X 10 

1,013 X 10^ 

k 

2.026 X 10 

1.013 X 10^ 

4 

1.013 X 10 


1.013 X 10^ 

2.026 X 10^ 

300 

300 

300 

300 

266 

266 

266 

266 

310 

310 

310 

310 


.830128 

.915968 


.514679 


.849515 

.982169 


.751^904 

.904845 





Table 2 

DETAILS OF ATMOSPHERIC MODELS 


I. R. -Active Constituent Concentrations (ppm) 


2 

Pressure (N/m ) 

Temperature (K) 

9.ll*2 X 10“ 

282 . 02 

7.136 

269.11 

5.169 

253.34 

3.176 

231.26 

9.555 X 10^ 

216.97 

I4. 

9.048 X 10 

291.30 

7.126 

280.00 

5.224 

264.80 

3.279 

242.80 

1.164 

218.60 


N^O CO SO^ CH, 


330 , 625 ^+. 

330 . 3 ^ 08 , 

330 . 1335 . 

330 . 21 k, 

330 . 15 . 


330 . li+ 310 . 
330 . 6529 . 

330 . 21U3. 

330. 50I+.3 

330. 17.84 
















APPENDIX A 


Deviations in inferred concentration as function of atmospheric model, species, 
cell partial pressure, and pollutant concentration, and deviations in 


Model 

Species 

Cell Partial 
Pressure 

(N/m^) 

Standard 

Pollutant 

Concent rat i on ( ppm ) 



A T. 
A 

Percent 



1 

H 

-.1 

-.01 

+ .01 

■a 

+1 

1 

CO 

, ^ h 

3. 5^+6 X 10 

.05 

- 80.0 

-8.0 

-1.0 

3.0 


50.0 




.1 

-TO.O 

-6.0 

-.5 

5.0 


1+2.0 




.2 

-66.0 

-5.0 

-1.3 

1.0 


1+0.0 




.3 

-66.7 

-5o3 

-.7 

.7 

Bi 

38.7 



1.013 X 10^ 

.05 

-152.0 


-2.0 

2.0 

D 

70.0 




.1 

-lUO.O 


-2.0 

.5 


60,0 




.2 

-177.0 

-10.3 

-.5 

1.0 

7.0 

56,0 




.3 

-217.0 

-9.3 

-.3 

2.0 

9.7 

57.0 



1.013 X 10^ 

.005 

a 

-70.0 

- 8.0 

i+.o 

68.0 





.0075 

* 

-U.O . 

- 5.3 

2.7 

U2.7 





.01 


- 28.0 

- 2.0 

k.O 

35.0 

* 




.015 

* 

-25.3 

- 2.7 

^ . 

2.0 

26.7 

* 












CH,. 

1.013 X 10^ 

,5 








4 



- 28.0 

- 5.0 

- 2.0 

2.0 


l8.0 ■ 




1.0 











- 27.0 

- 3.0 

- 1.0 

1.0 

3.0 

20.0 

























APPENDIX A (Cent. ) 


Model 

Species 

Cell Partial 
Pressure 

(N/m^) 

Standard 

Pollutant 

Concentration (ppm) 

At. Percent 
A 

-1 

-.1 

-.01 

+ .01 



rntm 

+1 

1 

CH. 

1.013 X 10^ 

1.5 

- 26.7 

- 2.7 

- .7 

■■ 


20.0 

(Cont, ) 

(cont *d) 












2.0 

- 32.0 

- 2.0 

- .5 


2.0 

20.5 




.5 

- 76.0 

- 8.0 

" ' 

-i+.O 

2.0 

ii.O 

1+2.0 




1.0 

-102.0 

- 7.0 

-2.0 

1.0 

k.o 

Ul.O 




1.5 

-133.0 

- 7.3 

- .7 

1.3 

6.1 

U.O 




2.0 

-110.0 

- 5.0 

- 1.0 

1.0 

8.p_ 

1+8.0 


NH3 

2.026 X 10^ 

.005 

- 88.0 

' - 8.0 

- 2,0 

' V . V ' 

~’67o 

- 

78.0 




.01 

- 52.0 

- k ,0 

- i.o 

1.0 

5.0 

1+5.0 




.02 

- 31.5 

- 2.5 

- .5 

•5 


28.0 


1.013 X 10' 


.025 

.005 

,01 

.02 


- 2Q,k 


- 2.8 

- .1+ 

.1+ 

-1+0.0 

- l+.o ' 

2.0 

-21+.0 

- 2.0 

3.0 

-11+.5 

- 1.5 

2.0 


2.h 

38.0 

26.0 


2h,h 


025 


-13.2 


.8 


1.6 


19.5 

17.6 










































APPENDIX A (Cont. ) 




Cell Partial 

Standard 

A T, Percent 
A 

Model 

Species 

(N/m^) 

Pollutant 

Concentration (ppm) 

-1 

■B 

-.01 

+ .01 

m' 

+1 

2 

CO 

3.5^6 X 10^ 

.05 

-l+O.O 

-U.o 


2.0 

u.o 

28.0 




.1 

-33.0 

-3.0 


2.0 

U.O 

2I1.0 




.2 

-32.0 

-2.0 

-.5 

1.0 

3.0 

2lt.5 




.3 

-31.7 

-2.0 

-.7 

1.3 

3.0 

18.0 


CO 

1.013 X 10^ 

.05 

-68.0 

-U.O 


2.0 

12.0 

Uii,0 




.1 

-6U.0 

-7.0 


1.0 

3.0 

38.0 




.2 

-70.0 

-8.0 


.5 

1.0 

36.0 




.3 

-90.0 

-3.3 

.3 

1.7 

7.0 

uo.o 


•SOg 

1.013 X 10^ 

.005 

-292.0 

-30.0 

- U.O 

2.0 

2I+.0 

* 

' 



.0075 

-19^.7 

-18.7 

- 2.7 

1.3 

17.3 





.01 

-ill 8.0 

-15.0 

- 2.0 

1.0 

13.0 

* 




.015 

-100.7 

-10.7 

- 1.3 

.7 

8.J 



CH^ 

1.013 X 10^ 

.5 

- 18.0 

- U.O 

- 2.0 

2.0 

h.O 

12.0 




1.0 

- 18.0 

- 3.0 

- 1.0 

1.0 

2.0 

12.0 




1.5 

- 16.0 

- 2.0 
1 

- .7 

.7 


13.3 






















APPENDIX A (Concl. ) 


Model 

Species 

Cell Partial 
Pressure 

Standard 

Pollutant 

Concentration (ppm) 

A Percent 

A 

-1 

Bl 

• 

o 

H 

+ .01 

D 

+1 

2 

CHi^ 

1.013 X 10^ 

2,0 


WB 

Bn 

H9I 

■■ 

999 

(Cont, ) 

( Cont . ) 

2.026 X 10^ 

.5 

- i+i+.o 

- 6.0 

- 2.0 

2.0 

h.O 

30.0, 




1.0 . 

- 53.0 

< k.o 

- 1.0„, 

...1.0 

h.O 

28.0 




1.5 

- 62.7 

- l+.T 

- 1.3 


h,l 

31.3 




2,0 

- 78.5 


- 1.0 


IBS 

324.0 



ii 

2.026 X 10 

.005 

- 3i^.0 

- 6.0 


g 

k.o 

32.0 




.01 

0 

1 — 1 
CVI 

1 

- 3.0 




19.0 




.02 

- 13.0 

- 1.5 


H 

Hi 

12.0 




.025 

- 11.6 

- 1.2 


B 

.8 

10.8 



1.013 X 10^ 

.005 


-l8.0 

- 2.0 

2.0 

lU.O 

* 




.01 

* 

- 9.0 

- 1.0 

2.0 

10. 0 

* 




.02 

» 

- 6.5 

- 1.0 

1.0 

5.5 

* 




.025 

* 

- 5.6 

- .k 

l.., M 

■ .k 

6.0 



^For these vaJ.ues of the parameters, the response functions of the instrument remained 
below the radiance levels used to infer concentrations. 
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APPENDIX B 


c(v) 


C'(v) 


E(v) 


N 


AN 


V 

AV 


Tq(v) 

V 

Av 


SYMBOLS 

transfer function of the total optical path through the gas 
cell at wavenumber v 

transfer function of the total optical path through the vacuum 
cell at wavenumber v 

radiance incident on the instrument at wavenumber v 

2 —1 
(watts /cm • sr cm ) 

average of radiances passing through the vacuum cell and gas 

o 

cell, integrated over wavenumber (watts /cm sr) 
difference between radiances passing through vacuum cell and 

p 

gas cell integrated over wavenumber (watts/cm sr) 
instrument output signal corresponding to N (volts) 
instrument output signal corresponding to AN (volts) 
transmissivity of the vacuum cell aperture, optical balance 
parameter (see fig. l) 

transmittance of the gas cell at wavenumber v 
wavenximber (cm”^) 

bandpass of the filter (cm~^) (see fig. 1) 


li+ 




Figure !• - Schematic of typical Gas Filter Correlation 
Radiometer* 
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I Figure 2* ■— Calculated response of AH signal in 

I CO chaimel to CO in U. S. Standard Atmosphere* 

”2 L 1962 9 for indicated partial pressures of CO in 

Instrument gas cell* 
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Figure 3. — Calculated response of aN signal in 
CHj^ ohsnnel to CHj^ in U. S. Standard Atmosphere, 

1962, for Indicated partial pressures of CH|^ in 

instrument gas cell. 
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Figure k, — Calculated response of AN signal in 
SOg channel to SO^ in U. S. Standard Atmosphere, 

1962, for indicated parti€U. pressure of SO^ in 

instrument gas cell. 
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Figure -- Calculated response of AH signal in 
channel to in S* Standard Atmosphere, 

1962, for indicated partial pressures of NH^ in 
instrument gas cell* 



PERCENT DEVIATION IN INFERRED CONCENTRATION 
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FiguT6 6 • ~ Positiv"© and. nsga"tiv6 
deviations in inferred CO^ 
concentration versus positive and 
negative deviations in for 

cells 1 and 2, in the U. S. Standard 
Atmosphere, 1962. 

(a) -Positive deviations, cell 1. 
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Figure 8.- Positive and negative 
deviations in inferred SO 2 
concentration versus positive and 
negative deviation in foi* 

cell 1, in the U. S. Standard 
Atmosphere , 19^2 . 

(a) Positive deviations, cell 1. 
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